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AERONAUTICAL SYMBOLS 
I. FUNDAMENTAL AND DERIVED UNITS 



Length. 
Time... 
Force-- 



Power- 
Speed _ 



Symbol 



I 
t 

F 



Metric 



Unit 



meter 

second 

weight of one kilogram . 



kg/m/sec. 
rkm/hr... 
I m/sec 



Symbol 



m 

sec 
kg 



English 



Unit 



foot (or mile) 

second (or hour) 

weight of one pound 

horsepower 

mi./hr 

ft./sec 



Symbol 



ft. (or mi.) 
sec. (or hr.) 
lb. 



HP. 

M. P. H. 

f. p. 8. 



TP, Weight, = mg 

g, Standard acceleration or gravity = 9.80665 
m/sec. 3 = 32.1740 ft./sec. 3 
W 

m, Mass, = — 
9 

p, Density (mass per unit volume). 
Standard density of dry air, 0.12497 (kg-m -4 

sec. 3 ) at 15° C and 760 mm = 0.002378 (lb.~ 

ft." 4 sec. 2 ). 

Specific weight of "standard" air, 1.2255 



2. GENERAL SYMBOLS, ETC 

ink 2 , Moment of inertia (indicate axis of the 
radius of gyration, h, by proper sub- 
script) . 



s, 


Area. 


s m 


Wing area, etc. 




Gap. 




Span. 




Chord length. 



kg/m 3 = 0.07651 lb./ft. 3 



Vj True air speed. 



6/c, Aspect ratio. 

Distance from c. g. to elevator hinge. 
Coefficient of viscosity. 



AERODYNAMICAL SYMBOLS 

7, Dihedral angle. 



q, Dynamic (or impact) pressure = ^ pF* 

L, Lift, absolute coefficient ^ = ~g 

D, Drag, absolute coefficient Cp^^ 

C y Cross - wind force, absolute coefficient 

Cc ~qS 

R, Resultant force. (Note that these coeffi- 
cients are twice as large as the old co- 
efficients L c , D c .) 

i w Angle of setting of wings (relative to thrust 
line). 

it, Angle of stabilizer setting with reference to 
to thrust line. 



VI Reynolds Number, where I is a linear 
p n 9 dimension. 

e. g., for a model airfoil 3 in. chord, 100 
mi./hr. normal pressure, 0° C: 255,000 
and at 15° C, 230,000; 
or for a model of 10 cm chord 40 m/sec, 
corresponding numbers are 299,000 
and 270,000. 
Center of pressure coefficient (ratio of 
distance of C. P. from leading edge to 
chord length). 
Angle of stabilizer setting with reference 

to lower wing, = (i t — i w ) . 
Angle of attack. 
Angle of downwash. 
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FLO'V AND DRAG FORMULAS FOR SIMPLE QUADRICS. 

ERRATA 

Paue 13, line 1, should rea.d: 

In the flow field q, P are constant rhere q^ 2 + q n 2 = 

constant . 

Pa^re 16, line 13, syrnhols should read: 

q x ,q y 

Page 33, line 3, Table III, should read: 
(For all shapes (l + Vy, 

q t = . ( 1 + i%) V sin g 

q n = - ( 1 - %5 V cos e ) 

Pan;e 23, Tatle III, delete line under table and substitute: 

*e is the angle between b ! and any normal to the coniocal 
surface . 

Pmgtt 23, Table HI: 

cos a) should be omitted from formulas. 
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FLOW AND DRAG FORMULAS FOR SIMPLE QUADRICS 

By A. F. Zahm 



PREFACE 

In this text are given the pressure distribution and resistance found by theory and experi- 
ment for simple quadrics fixed in an infinite uniform stream of practically incompressible fluid. 
The experimental values pertain to air and some liquids, especially water; the theoretical refer 
sometimes to perfect, again to viscid fluids. For the cases treated the concordance of theory 
and measurement is so close as to make a resume of results desirable. Incidentally formulas 
for the velocity at all points of the flow field are given, some being new forms for ready use 
derived in a previous paper and given in Tables I, III. A summary is given on page 22. 

The computations and diagrams were made by Mr. F. A. Louden. The present text is a 
slightly revised and extended form of Report No. 312, prepared by the writer for the Bureau of 
Aeronautics in June, 1926, and by it released for publication by the National Advisory Com- 
mittee for Aeronautics. A list of symbols follows the text. 

PRESSURE AND PRESSURE DRAG 

We assume the fluid, of constant density and unaffected by weight or viscosity, to have in 
all the distant field a uniform velocity q 0 parallel to x; in the near field the resultant velocity q. 
If now the distant pressure is everywhere p 0 , and the pressure at any point in the disturbed flow 
is po + V, the superstream pressure p is given by Bernoulli's formula, 

VlVn = l~tflq.o\ (1) 

where p n = pq 0 2 /2, called the "stop" or "stagnation " or "nose" pressure. 

At any surface element the superpressure exerts the drag JT p dy dz, whose integral over 
any zone 1 of the surface is the zonal pressure drag, 

D = fpdydz. (2) 

Values of p } D are here derived for various solid forms and compared with those found by 
experiment. 

PRESSURE MEASUREMENTS 

The measured pressures here plotted were obtained from some tests by Mr. R. H. Smith 
and myself in the United States Navy 8-foot wind tunnel at 40 miles an hour. Very accurate 
models of brass, or faced with brass, had numerous fine perforations, one at the nose, others 
further aft, which could be joined in pairs to a manometer through fine tubing. Thus the 
pressure difference between the nose and each after hole could be observed for any wind speed. 
Then a fine tube with closed tip and static side holes was held along stream at many points 
abreast of the model, to show the difference of pressure there and at the nose. Next the tube 
was thrust right through the model, to find the static pressure before and behind it. The 
method is too well known to require further description. 

THE SPHERE 

Assume as the fixed body a sphere, of radius a, in a uniform stream of inviscid liquid, as 
shown in Table I. Then by that table the flow speeds at points on the axis x, y and on the 
surface are 

q x = (1 -a?/^, q v = (1 + a z /2y>)a 0 , 2< = 1.52, sin 0, (3) 

where 6 is the polar angle. Figure 1 shows plots of these equations. 



i a zone is a part of the surface bounded by two planes normal to qo. Usually one plane is assumed tangent to the surface at its upstream end. 

3 
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To graph p/p n in Figure 1, we subtract from the line y = 1, first q x 2 /q 0 2 to show the pressure 
along x; then q_t\q_o to portray the surface pressure. A similar procedure gives the superpressure 
in the equatorial plane. 

The little circles show the actual superpressures found with a 2-inch brass sphere in a tunnel 
wind at 40 miles an hour. These agree well with the computed pressures except where or 
near where the flow is naturally turbulent. 

By (3) and (1), on the sphere's surface ^/p n = 1-2.25 sin 2 0; hence the zonal pressure drag 
fp.2wydy is 

7ra 2 sin 2 0(l- |sin 2 0)p n , (4) 



D 



for a nose cap whose polar angle is 0. With increase of 0, as in Figure 2, D/p n increases to a 
maximum .698- a 2 for 0 = 41° -50' and p = 0; then decreases to zero for 0 = 70° -37'; then to its 
minimum — .3927 a 2 for 9=t/2; then continues aft of the equator symmetrical with its fore 
part. Thus the drag is decidedly upstream on the front half and equally downstream on the 




r 

3 * 5 

Length in inches 



Wind 



Fig. 1.— Velocity and pressure along axes and over surface of sphere; graphs indicate theoretical values; 
circles indicate pressures measured at 40 miles per hour in 8-foot wind tunnel, United States Navy 

rear half, having zero resultant. The little crosses, giving D/p n for the measured pressures, 
show that the total pressure drag in air is downstream, and fairly large for a body so blunt as 
the sphere. 

Figure 3 depicts the whole-drag coefficient 2 G D = 2D/ir p a 2 q 0 2 , of a sphere, for the manifold 
experimental conditions specified in the diagram, plotted against Reynolds Number R = 2 qoa/v, 
v being the kinematic viscosity. For 0.2 < R < 200000, the data lie close to the line. 



Cz> = 28£-- 85 + .48, 

an empirical formula devised by the writer as an approximation. 
For .5<i2<2 (5) fairly merges with Oseen's formula 

C D = 24#~ 1 +4.5, 



(5) 



(6) 



and for R < .2 Stokes' equation C D = 24/# is exactly verified. Both these formulas are theoretical. 
Stokes treated only viscous resistance at small scale; Oseen added to Stokes' drag coefficient, 
24/#, the term 4.5 due to inertia. 

3 From the drag D= Cp n S, where S is the model's frontal area, one derives the drag coefficient Cd=D/P*S. 
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Over an important R range Figure 3 shows C D = .5, giving as the sphere's whole drag 

D = .5p n S, (7) 

where S = w a 2 is the frontal area. That is, the sphere's drag equals half its nose pressure times 
its frontal area, For R<.2 Stokes' value, D = 6t ;u a q n , has been exactly verified experi- 
mentally, as is well known. 



y 




\ =1.0! 
\ ' 

N ' 

—1.5 

Fig. 2.— Pressure and pressure-drag on sphere. Graphs indicate 
theoretical values; circles indicate pressure p/p n measured at 
40 miles per hour; crosses indicate pressure-drag -D/pn, com- 
puted from measured pressure 

THE ROUND CYLINDER 

Next assume an endless circular cylinder, of radius a, fixed transverse to the stream, as 
indicated in Table I. By that table the flow speed at points on the axes x, y and on the surface is 

g x = (1 - a 2 /x 2 ) q 0 , (1 + a 2 fy 2 ) g 0 , it = 2 So sin d y (8) 

where 6 is the polar angle. Plots of (8) are shown in Figure 4. 

Graphs of p/p n , made as explained for the sphere, are also given there, together with experi- 
mental values, marked by small circles, for an endless 2-inch cylinder in a tunnel wind at 40 
miles an hour. The agreement is good for points well within the smooth-flow region. 

On the surface p/p n =l-4: sin 2 d. The integral 2 f pdy gives, per unit length of 
cylinder, the zonal pressure-drag formula, 

Q 

D/p n = 2 a sin 0-^ a sin 3 0. (9) 



SYMBOLS:- 

B ■ Drag of sphere 
d = Diam. - 

S = Areo » - nd*/4 

V = Relative speed of sphere & fluid 
p = Density of fluid 

V m Kinematic viscosity 



. h. S.Allen 
> ' Phil. Mag, 1900 




Eiffel. Sphere in air. d = 16.2 cm 
"A/ouvel/es Recherches" etc ., IS 14. 
0. W. Silvey. Mercury drops in castor 

oil. d'.Om to .0756 cm. Rhys., Rev, 19/6. 
H. D. Arnold. Rose metal spheres /n 

co/zo oil. d = .OI3to. 141 cm. Phil. Mag., 191 L 
Liebster & Schiller. Steel spheres in 

glycerin, sugar solutions & water. 

d=.l to Jem. Rhys. Zeit, 1924. 
Air bubbles in woier. > 

d=.0094 to .06/ cm 
Air bubbles in aniiine 

d = .007 to .110 cm 
Paraffin spheres in aniline 

d = .069 to. 3/6 cm 
Amber spheres in water 

d = .l/4 to .346 cm 
Steel spheres in water 

d'. 3/6 to .792 cm 
d ' * .8 cm ~\ Steel spheres ^ 
d = 1.8 cm j in air „ ,J- , . 
d = 6.21 cm < C.Wieselsberger 

d = 9.98 cm I, Hollow copper > y ,qpp 

d = /4.2 cm [spheres in air 

d = 28.25 cm) ) 

A. F. Zohm. Polished pine sphere in air 

( Without correction for air compression) 

d = 10.12 cm (7620 < V< 27432 cm/sec.) 

RhiL.Mag., 1901. 
Atmospheric Tunnel ""j Maple spheres in 

d= 15 & 20 cm I air. 0. L. Bacon 

Variable Density Tunnel \ and t~. G.Reid 

d=20 cm J N.A.C.A. Rep't.l85 
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12 5 10 
Reynolds Number, R = Vd/u 

Fig. 3.— Drag coefficient for a sphere in steady translation through a viscous tluid 




Fig. 4 —Velocity and pressures along axes and over surface of endless cylinder; graphs indicate theoretical values; 
circles indicate pressures measured at 40 miles per hour in 8-foot wind tunnel, United States Navy 



Wind 



O Infinitely long wires in air, 
d - . 0056 cm io 3. 1 75 cm, 
Aeronautics Staff of N.P.L., 
British R.&M. No. 102 
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C D = SnlR/2. 002 - log e R) 
/ L omb, Phil. Mag., 1911 



C D = 9.4R- d + /.2 



= .005 crrP\ 
.01 
.03 
.10 
.30 
.19 
4.2 

8.0 

30. 



Infinitely long ^ 
cylinders in air. 
* C.Wieselsberger 
Rhys. Zeit, 1921 



D = Drag of cylinder 

per diameter length 
d = Diameter of cylinder 
S = frontal area per 

diameter-length * a 
V= Relative speed of 

cylinder and fluid 
p = Density of fluid 
V = Kinematic viscosity 



IO 4 



IO 3 



-3.0 ±' 



Fig. 5.— Pressure and pressure-drag on endless 
cylinder. Graphs indicate theoretical values; 
circles indicate pressure vlPn measured at 40 
miles per hour; crosses indicate pressure-drag 
DJp n computed from measured pressure 



10 20 50 10 2 IO 3 

Reynolds Number, R = Vdlv 

Fig. 6.-Drag coefficientf or an endless cylinder in steady translation through a viscous fluid 
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This is 0, 2a/3 (max.),0, -2a/3, for 0 = 0°, 30°, 60°, 90°; and is symmetrical about the equatorial 
plane z = 0. In Figure 5, the little crosses give D/p n for the measured pressures, and show total 
D/p n = 2.33a. 

Figure 6 delineates the drag coefficient G D plotted against R = 2 aq 0 /v, from Wieselsberger's 
(Reference 1) wind tunnel tests of nine endless cylinders held transverse to the steady flow. 
The faired line is the graph of 

(7 D = 9.4fl-- 8 +1.2, (10) 

an empirical equation devised by the present writer. 

For very low values of R, Lamb derives the formula 



Ud (2.002 -lo&fl)fi' 
whose graph in Figure 6 nearly merges with (10) at R= .3. 

For 1 5000 <R< 200000, Figure 6 gives Co ~ 1.2; hence the drag per unit frontal area is 

D~1.2p nl (12) 
which is 2.4 times that for the sphere, given by (7). 



THE ELLIPTIC CYLINDER 

An endless elliptic cylinder held transverse to the stream, us shown in Table I, gives for 
points on x, y and on its surface, 

q x -(l-n)g 0f £„« q t = (1 + b/a)q 0 sin 0, (13) 

where m, n are as in Table I. Amidships q t = (1 + b/a)q 0 = 2q 0 for a = b, as given by (8). Graphs 
of (13) are given in Figure 7. 

To find a', V for plotting (13), assume a' and with it as radius strike about the focus an 
arc cutting y. The cutting point is distant V from the origin. Otherwise, V = ^a n - c 2 , where 
c 2 = a 2 — & 2 = const. 

With a/b = 4 one plots p/p n in Figure 7, as explained for the sphere. The circles give the 
experimental p/p n for an endless 2-inch by 8-inch strut, at zero pitch and yaw, in a tunnel wind at 
40 miles an hour. The theoretical and measured pressures agree nicely for all points before, 
abreast, and well behind the cylinder. 

Again, sin 2 ^ = aV/(& 4 + cV), if c 2 = a 2 -6 2 . Hence on the model 

p/Pn=l-qtVao 2 =l- { j^f' (14) 

y 

This gives the zonal pressure drag, D = 2f pdy, per unit length of cylinder, or 



^^2,-2(a +6 )/^ 2 = -4 6 a + ^ + 2 & ^ 2 tan-'|, (15) 
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whose graph, for a/6 =4, appears in Figure 8. 
p = 0, then falls to its minimum amidships. 



It rises from 0 at the nose to its maximum where 




Fir,. 7. 



t — r~ 

6 p 10 l- 

W/nd - Length in inches 

-Velocity and pressure along axes and over surface of endless elliptic cylinder. Graphs indicate theoretical values; 
circles indicate pressure measured at 40 miles per hour in 8-foot wind tunnel, United States Navy 




Fig. 8.— Pressure and pressure-drag on endless elliptic cylinder. Graphs indicate theoretical values; circles 
indicate pressure p/p» measured at 40 miles per hour; crosses indicate pressure-drag D/p n computed from 
measured pressure 

Whatever the value of a/b, the whole pressure on the front half is negative or upstream, 
as for the sphere and round cylinder, and is balanced by the rear drag. For b fixed it decreases 
indefinitely with b/a. 
10420—27 2 
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The crosses marking actual values of D/p n found in said test show a downstream resultant 
D. In fact, it is one-third the whole measured drag of pressure plus friction, or one-half the 
friction drag. 

For the cylinder held broadside on, b>a and a 2 -b 2 = -c 2 , hence changing c 2 to -c 2 under 
the integral sign of (15), we find 



(16) 



where now c 2 = b 2 -a 2 . With b fixed, the upstream pressure drag on the front half increases with 
b/a, becoming infinite for a thin flat plate. It is balanced by a symmetrical drag back of the 
plate. 

Such infinite forces imply infinite pressure change at the edges where, as is well known, the 
velocity can be q = -J2p r /p = oo , in a perfect liquid whose reservoir pressure is p r = co . Otherwise 
viewed, the pressure is p r at the plate's center (front and back) and decreases indefinitely toward 
the edges, thus exerting an infinite upstream push on the back and a symmetrical downstream 
push on the front. In natural fluids no such condition can exist. 

THE PROLATE SPHEROID 

A prolate spheroid, fixed as in Table I, gives for points on x, y and the solid surface, respec 
tively, the flow speeds 



&=(l-7&)&, 



(l + m)q 0 , 



(14- & a ) q 0 sin 0, 



-to 




(16) 



Wind 



1 r t 
6 8 10 

^Length in inches 



Fig. 9.— Velocity and pressure along axes and over surface of prolate spheroid. Graphs indicate theoretical values; circles 
indicate pressures measured at 40 miles per hour in 8-foot wind tunnel, United States Navy; dots give pressures found 
with an equal model in British test, R. and M. No. 600, British Advisory Committee for Aeronautics 




Fig. 10.— Pressure and pressure-drag on prolate spheroid. Graphs indicate theoretical values; dots indicate 
measured pressure p/Pn from Figure 9; crosses indicate pressure-drag Dfp n computed from me;isured pressure 
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where Tc a is to be taken from Table II. Graphs of (16) are given in Figure 9, for a model having 
a/b = 4, viz., ~k a = 0.082. 

For this surface p/p n plots as in Figure 10. For a 2 by 8 inch brass model values of p/p n 
are shown by circles for a test at 40 miles an hour in the United States Navy tunnel; by dots for 
a like test in a British tunnel. (Reference 2.) 

By (16), for points on the surface p/p n = 1 -q t Vio 2 = 1 ~ (1 sin 2 6. From this, since 

sin 2 d = a 2 y 2 / (b A + c 2 y 2 ), the zonal pressure drag / p. 2 tt y dy is found. Thus 



Dlp n = iry 2 



^(l + lc a ) 2 y 2 + c , 



^*&*/i_±_i vji h' + c 2 y 2 

(l+fc a ) 2 l0g g — . 



(17) 



Starting from y = 0, D/p n increases to its maximum when p = 0, or sin 0=l/(l + fc a ); then 
diminishes to its minhnum for y = b. Figure 10 gives the theoretical and empirical graphs of 
Djp n for a/b = 4. 

For b fixed the upstream drag on the front half decreases indefinitely with b/a, becoming 
zero for infinite elongation. 

OBLATE SPHEROID 

The flow velocity about an oblate spheroid with its polar axis along stream is given by 
formulas in Table I, and plotted in Figure 11, together with computed values of pjp n . No 
determinations of p or D were made for an actual flow. The formula for D/p n is like (17), 
except that c 2 = b 2 -a 2 , and Jc a is larger for the oblate spheroid, as seen in Table II. For b fixed 
the upstream drag on the front half increases indefinitely with b/a. 




-e 

Wind 

Fkj. 1 1 — Theoretical velocity and pressure along X axis of oblate spheroid. Diameter/thickness =4 

CIRCULAR DISK 

The theoretical flow speeds and superpressures for points on the axis of a circular disk fixed 
normal to a uniform stream of in viscid liquid are plotted in Figure 12, without comparative 
data from a test. One notes that the formulas are those for an oblate spheroid with eccen- 
tricity 6=1. . 

For 1500<g 0 a/* <500000, Wieselsberger (Reference 3) finds for the air drag ol a thin 

normal disk, of area S, 

D=l.l p n S, (18) 
or 2.2 times that for a sphere. For Ofy/v extremely small, theory gives 

D = 5.1 tt ix a q 0 , (19) 
; weU known. Test data for a complete graph, including these extremes, are not yet available. 
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Fig. 12— Theoretical pressure and velocity along axis of disk 
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REGIONS OF EQUAL SPEED 
In the flow field q, p are constant where <Zt 2 + & 2 = constant, viz. where 
2 2 /g 0 2 =(l+ra) 2 sin 2 0+(l-n) 2 cos 2 0 = const. 
In particular for the region q = q 0 , this becomes 



tan 2 0 = - 



n 2 — n a' 



m2 + m &' 4 



tan 2 0* 



(20) 



(21) 
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Fig. 13 —Inertia coefficient vs. elongation. Plotted from Table II 



O 



3* 

.Oj 



to 



which applies to all the quadrics in Table I. Clearly tan 0 = 0 for n = 2; ta,n 2 $ = n/m for 
m, n = 0, viz. for all distant points of (21). For these points the normal to any confocal ellipse 
lies along the radius vector and asymptote of (21), as seen in Figures 14 to 17. 




Fig. 14.— Lines of steady flow, lines of constant speed and pressure, for 
infinite frictionless liquid streaming past a sphere 

For the sphere n = 2m = a?/r z ; hence (21) becomes 

tan 2 0= 2^^5-3 
where r = a' = ~Jx 2 + y 2 . The form of this is depicted in Figure 14. 



(22) 



* tan p=y/x is the slope of a radial line through the point (x, y) where (21) cuts a confocal curve a'V, of Table I. Knowing a', &', 0, to locate 
(j, y) draw across the radial line an arc of a'b' by sliding along the x, y axes a straightedge subdivided as in the ellipsograph. The operation is rapid 
and easy. 
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For a round cylinder n = m, = a 2 /r 2 ; hence 

tan 2 0= ^ 9-. — o or, 2r 2 = a 2 sec 26. 
2r 2 + a 2 1 ' 

which is the section of a hyperbolic cylinder, as in Figure 15. 

V 




Fig. 15 —Lines of steady flow, lines of constant speed and pressure, for infinite frictionless liquid streaming across endless 

round cylinder 




Fig. 16.— Lines of steady flow, lines of constant speed and pressure, for infinite frictionless liquid 
streaming across endless elliptic cylinder 
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A plot of (21) for an elliptic cylinder, fixed as shown in Table I, is given in Figure 16; for a 
prolate spheroid in Figure 17. 

Besides the region (21) , having q = q 0 , it is useful to know the limit of perceptible disturbance 
say where q 2 /q 0 2 =l ± .01. This in (20) gives 



(1 + m) 2 sin 2 6+ (1-n) 2 cos 2 0=1 ±.01, 
which applies to all the quadrics here studied. Hence 

n 2-n 0.01 



tan 2 6 = 



m 2 + m ra(2 + m)cos 2 0 



(24) 



(25) 



A graph of (25) for a round cylinder is shown in Figure 15. Like plots for the other quadrics 
are readily made. 




Fig. 17— Lines of steady flow, lines of constant speed and pressure, for infinite friction- 
less liquid streaming past a prolate spheroid. Full-line curve q=q« refers to stream 
parallel to x; dotted curve q=q 0 refers to stream inclined 10° to x. 

If in (20) a series of constants be written for the right member, the graphs compose a family 
of lines of equal velocity and pressure, covering the entire flow field. Rotating Figures 14, 17 
about x gives surfaces of q = q 0 - 

COMPARISON OF SPEEDS 

Before all the fixed models the flow speed is q 0 at a great distance and 0 at the nose; abreast 
of them it is q 0 at a distance, and (H-& a )g 0 amidships. 

The flux of q — q 0 through the equatorial plane obviously must equal q 0 S where S is the 
body's frontal area. Hence two bodies having equal equators have the same flux q 0 S, and the 
same average superspeed or average q-q 0 . But the longer one has the lesser midship speed; 
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hence its outboard speed wanes less rapidly with distance along y. A like relation obtains 
along x from the nose forward. These relations are shown in the velocity graphs of Figures 
18 and 19. A figure similar to 18, including many models, is given in Reference 4. 




Fig. 18.— Superposed graphs of How speed abreast of endless round and elliptic 
cylinders of same thickness fixed transverse to an infinite stream of inviscid 
liquid. At great distance flow speed is q 0 




1 10. 19, — Superposed graphs of axial flow speed before three endless cylinders 1, 2, and 3 (3 escalating 
2) , each fixed transverse to an infinite stream of inviscid liquid. At great distance flow speed is q 0 

COMPARISON OF PRESSURES 



The foregoing speed relations determine those of the pressures. The nose pressures all 
are p n = pq 0 2 /2; the midship ones are p=p n — (1 +Jc a ) 2 p n ' The drag on the front half of the 
model is upstream for all the quadrics here treated; it increases with the flatness, as one proves 
by (15), (17), and is infinite for the normal disk and rectangle. 

APPLICATION OF FORMULAS 

The ready equations here given, aside from their academic interest in predicting natural 
phenomena from pure theory, are found useful in the design of air and water craft. The formula 
for nose pressure long has been used. That for pressure on a prolate spheroid, of form suitable 
for an airship bow, is so trustworthy as to obviate the need for pressure-distribution measure- 
ments on such shapes. The same may be said of the fore part of well-formed torpedoes deeply 
submerged. The computations for stiffening the fore part of airship hulls can be safely based 
on theoretical estimates of the local pressures. The velocity change, well away from the 
model, especially forward of the equatorial plane, can be found more accurately by theory 
than by experiment. The equation (21) of undisturbed speed shows where to place anemometers 
to indicate, with least correction, the relative speed of model and general stream. 

REFERENCES 

1. Wieselsberger, C: Physicalische Zeitschrift, vol. 22. 1921. 

2. Jones, R., and Williams, D. H.: The distribution of pressure over the surface of airship model U. 721, 

together with a comparison with the pressure over a spheroid. Brit. Adv. Com. for Aeron. Reports and 
Memoranda No. 600. 1919. 

3. Wieselsberger, C: Physicalische Zeitschrift, vol. 23. 1922. 

4. Taylor, D. W.: Speed and Power of Ships, gives a figure similar to 18 but including more models. 1910. 
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SYMBOLS USED IN TEXT 



x, y Cartesian coordinates; also axes of same. 

r, /3 Polar coordinates. 

a Angle of attack of uniform stream. 

s Length of arc, increasing with ft. 

0 Inclination to x of normal to confocal 

curves in Table I. 

<p Velocity function. 

\J/ Stream function. 

q Resultant velocity at any point of fluid. 

q 0 Velocity of distant fluid (parallel to x 

axis). 

qcn Qr Velocity at points on x and y axes (parallel 

to x axis). 

q t Velocity along confocal surface or model 

surface. 

q n Velocity normal to confocal surface. 

p Density of fluid. 

n Viscosity. 



v Kinematic viscosity. 

p n Nose pressure =p qj-ll. 

p 0 Pressure in distant fluid. 

p Superstream pressure anywhere. 

D Zonal pressure drag=y* p dy dz. 

D Whole drag. 

S Frontal area of model. 

Cd Drag coefficient=Z) / p n $. 

R Reynolds number. 

a Radius of sphere, cylinder. 

a, b Semiaxes of ellipse. 

a', b' Semiaxes of confocal ellipse. 

e Eccentricity of ellipse. 

e' Eccentricity of confocal ellipse. 

c Focal distance = ae = a' e' = ja 2 — b 2 

k a Inertia factor (Table II). 

m, n, m 8 .~- Quantities defined in Tables I, II. 



TABLE I 

Flow functions for simple quadrics fixed in a uniform stream of speed q Q along x positive 



Symbol defini- 
tions 



See diagram A 
(fig. 20) 



Sec diagram B 
(fig. 20) 



See diagram C 
(fig. 20) 



Form of quadric 



Sphere 



Circular cylinder 



Elliptic cylinder 



Prolate spheroid 

e=~^a^-b 2 



Oblate spheroid 



Circular disk 
a = 0, e=l 



Value of functions at any confocal surfaces of semiaxes a' } b' 



Velocity function <p 



<P= — (1 + m) qoX, where 
a 3 

m= 2^ 



(p= — (l-\-m) qa, 



<P= — (l + m) qoX, 

_b a+b 
m ~a' a' + b' 



— (1 + ro) qoX, 
~~ 1+e 2e~ 

^g, jpi-iiqji 



<P= — {l + m) QoX, 
e'b' 



oi'n ' f> f 



m=— I —,— sin- l e' ) 



Stream function \p 



2 (l-n)q 0 y 2 , where 



a* 
a 3 



t=-(l-n) q a y, 



t=-(l-n) q 0 y, 

_b^ a+b 
n ~V a' + b' 



^ = — | (1-n) q 0 y 2 , 



log* 
log, 



1 + e' 2e' 
l-e' l-e'' 
1 + e 2e 2 



l-e l-e 2 



_1 
2 

e'a 



d-n)q 0 y 2 , 



ea 



■sin l e 



^=~2 (!— n ) Qoy 2 , 
2 



Component velocities 



Differentiation along 

arc s of either figure 

gives: 
d<p dx . N 
qt= dx ds =il + m)g » 

sin 0, valid for all the 

figures; 
„ tydy n s 

cos 0, for the cylin- 
ders; 
1 dxf/ dy 
dy 1= 
qo cos 0, for the axial 
surfaces; viz., sphere, 
spheroids, disk. 

For a', b' = a, b, Table 
II gives m 8 ; whence 

q t =(l + m 8 ) q 0 sin 0, 

as the flow velocity on a 
fixed quadric surface. 

q n = ±0 for disk, since 
n=l 

Remark — both q x q n 
can be derived from 
either <p or f'\ 

If Qt <7n = max. q t q n 
on a'b' , at any other 
point thereof 

q t = q t sin 0, q n =q n 
cos 0 



<p, \p, in elliptic coordinates, can be found in textbooks; e. g., §§ 71,105,108, Lamb's Hydrodynamics, 4th Ed. 
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TABLE II 



Inertia factors k a * for quadric surfaces in stead}- translation along axis a in Figure 20 



Elonga- 
tion E 


Elliptic cylinder, E=a/b 
a 


Prolate spheroid E — a/b 

14-c 
logc i — ; — 2c 
1 — e 

[ ° ge l-e 1-e* 


Oblate spheroid E h\n 

. en — n sin l e 

k(l ~~~~ e — E shr l e 


L 00 


1. 000 


0. 500 


0. 500 


1. 50 


. 667 


. 305 


. S03 


2. oo 


. 500 


. 209 


1. 118 


2. 50 


. 400 


. 157 


1. 428 


3. 00 


. 333 


. 121 


1. 742 


4. 00 


. 250 


. 082 


2. 379 


5. 00 


. 200 


. 059 


3. 000 


6. 00 


. 167 


. 045 


3. 642 


7. 00 


. 143 


. 036 


4. 279 


8. 00 


. 125 


. 029 


4. 915 


9. 00 


. Ill 


. 024 


5. 549 


10. 00 


. 100 


. 021 


6. 183 


00 


. 000 


. 000 


00 



* In this table k a = m t of Table I, viz, the value of m when a', b'=a, b. Lamb (R. and M. No. 623, Brit. Adv. Com. Aeron.) gives the numerical 
values in the t hird column above. For mot'on of elliptic cylinder along b axis inertia factor is kb=a/b. 



Diagram A y Diagram B y Diagram C y 




Kin. 2(1 



VELOCITY AND PRESSURE IN OBLIQUE FLOW ' 

PRINCIPLE OF VELOCITY COMPOSITION 

A stream q 0 oblique to a model can be resolved in chosen directions into component streams 
each having its individual velocity at any flow point, as in Figure 21. Combining the individuals 
gives their resultant, whence p is found. * 

VELOCITY FUNCTION 

Let a uniform infinite stream q 0 of inviscid liquid flowing past a fixed ellipsoid centered 
at the origin have components U, V, W along x, y, 2, taken parallel, respectively, to the semi- 
axes, a, b, c; then we find the velocity potential <p for q 0 as the sum of the potentials <p a , <p b , <p c 
for U, V, W. 

In the present notation textbooks prove, for any point (z, ?/, z) on the confocal ellipsoid 
a' V c' 

<p a =-(l + m a )Ux, (26) 

and give as constant for that surface 

* A f°° da ' V 1 f°° da/ * ,m 

m a = abc(l-abc ^y^) fl , 3 p ? (20 

/•oo 

the multiplier of j being constant for the model, and X = a' 2 -a 2 . Adding to (26) analogous 

values of <^ 6 , <p c gives 

0= - (!+ <m a ) Ux- (1 + m b ) Vy- (1 + m c ) Wz = - (1 + m)q 0 h, (28) 



2 This brief treatment of oblique flow was added by request after the preceding text was finished. 

* Simple formulas for this integral and the corresponding b, c ones, published by Greene, R. S. Ed. 1833, are given by Doctor Tuckerman in 
Report No. 210 of the National Advisory Committee for Aeronautics for 1925. Some ready values are listed in Tables III, IV. 
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where Ji is the distance of (x, y, z) from the plane <p = 0, and m a) m b , m c , m are generalized inertia 
coefficients of a' h' c' for the respective streams U, V 9 W, q 0 . For the model itself the inertia 
coefficients usually are written Jc a , Tc c , h. The direction cosines of h are 

£-\±28 5, Jf_l+2L» I, N= 1 ^ ^, (29) 
1 -f- ra g 0 1 -f ra <fr 1 + m q 0 v/ 

as appears on dividing (28) by (1 + m)q 0 , the resultant of (1 + rn a ) E7, (l + m b ) V, (1 4- m c ) W. 

EQUIPOTENTIALS AND STREAMLINES 

On a' V c r the plane sections <p = constant are equipotential ellipses parallel to the major 
section <p = 0, and dwindling fore and aft to mere points, which we call stream poles, where the 
plane (28) is tangent to a' V c' . If e is the angle between any normal to a' V c' and the polar 




normal, whose direction cosines are i, M, N, we call the line e = const, a line of stream latitude. 
Thus e is the colatitude or obliquity of a surface element of o! V c'. The line e = 90° is the stream 
equator. This latter marks the contact of a tangent cylinder parallel to the polar normal, viz, 
perpendicular to the plane (28), as in Figure 22. If Z, m, n are the direction cosines of any normal 
to a' V c' 

cos e = lL + mM+nN. (30) 

Since the streamlines all cut the equipotentials squarely, 3 the polar streamline must run 
continuously normal to the family of confocal ellipsoids a' V c' . Hence it forms the intersec- 
tion of a pair of confocal hyperboloids, and at infinity asymptotes a line parallel to q 0 through the 
origin. This straight line may be called the stream axis. Its equation is x: y: z - U:V:W. 



1 On the model, therefore, the streamlines are longitude lines, viz. orthogonals to the latitude lines. 
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COMPONENT VELOCITIES 

At any point of any confocal surface a' V c/ the streamline velocity q, perpendicular to the 
cquipotential ellipse there, has components q ni q tJ respectively, along the surface normal n 
and the tangent s in the plane of q and n. By (28) we have 

dp dh - . , Di v 
2 <~Mdi te2 ' 3me ' (31) 

where — d(p/dh = (1 4- m)g 0 =gi = max. ^, is the equatorial velocity. By (26) the inward normal 
velocity due to <p a is 

-J (l+m a )Ux=-l(l-n a )U, (32) 

O n 

n a being constant on a' b' c', as may be shown. Similarly, <p b , <p c contribute — m(l — n b ) V, — n 
(l — n c )W; hence the whole normal component is 

g B = -l(l-n a )U-m(l~n b )V-n(l-n c )W==q n cose, (33) 

where q n = [(l — n a ) 2 U 2 + (1 — n b ) 2 V 2 + (1 — n c ) 2 W 2 ]- 5 = ma,x. # n is the normal velocity at the stream 
poles. Some values of n a , n b are given in Tables I, III. One also may find (33) as the normal 
derivative of (28). 

We now state (28) : At any point of a' b' c' the velocity potential equals q t h, the equatorial 
speed times the distance from the plane of zero potential. Similarly (31) (33) state: At any 
point of a' b' c f the tangential speed (q t sin e) equals the equatorial speed times the sine of the 
obliquity; the normal speed (q n cos e) equals the polar speed times the cosine of the obliquity. 
This theorem applies to all the confocals, even at the model where q n = 0 4 . 

Incidentally the normal flux through a' V c' is fq n cos €• dS = fdSf, where is the pro- 
jection of S on the plane of </> = const, and equals the cross section of the tangent cylinder. 
The whole flux through a' b' c' is therefore zero, as should be. 

POLAR STREAMLINE 

Some of the Foregoing relations are portrayed in Figure 22 for a case of plane flow. Note- 
worthy is the polar streamline or hyperbola. Starting at infinity parallel to q 0 , the polar fila- 
ment runs with waning speed normally through the front poles of the successive confocal sur- 
faces; abuts on the model at its front pole, or stop point; spreads round to the rear pole; then 
accelerates downstream symmetric with its upstream part. Its equation q t = 0 = d<£/ds can be 
written from (28) 

q t = (l+m a )U sin 6— (l + m & )Fcos 6 = 0, or tan Q=\^^ft (34) 

This asymptotes the stream axis y/x = V/U; for at infinity m a , m b = 0, and tan d=V/U. Plane- 
flow values of m a , m b are given in Tables I, III. 

All the confocal poles are given by (34); those of the model are at the stops where 

, Q l + £a V a 2 y _ 
tane= iT¥ b U = b 2 x' W 

Thus on an elliptic cylinder they are where y/x = ¥/a 3 . V/U; on a thin lamina they are at 
x = ^ c C os a, as given in the footnote. Tables II, IV give values of Jc a , lc b . 

* An analogous theorem obtains also for any other uniform steady stream, say of heat or electricity, that has zero normal component at the bound- 
ary ellipsoid and zero concentration in the flow field. 

* To graph (34) we may use the known relations. f 

tan 0=l,* y x =p tana, (35) 

where tan a = V/U is the slope of q 0 or the asymptote to (34). Thus (34) becomes a7&'=(l+m&)/(l+m a ), which with the tabulated values of 
m a , mb, reduces to 

VJ- =1 (36) 

c«cos 2 « c 2 sin 2 « ' 

a hyperbola whose semiaxes are c cos a, c sin a, c being the focal distance. In this treatment x=a' cos a, y=V sin a, a being a fixed eccentric 
angle of the successive confocal elipses. 
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Each angle of attack has its own flow pattern; each its polar streamline given by (34). 
A close-graded family of confocal ellipses and hyperbolas therefore portrays all the poles and 
polar streamlines in the plane ab for all angles of attack. The family can be written 

y = V sin a. 



x = a cos a, 

Thus, giving a', W a set of fixed values, then a a set, we have the confocal families 



x 4. y = 1 



= 1. 



(38) 



(39) 



(40) 



the first being ellipses, the second hyperbolas like (36) below. 

Similarly, the locus q n = 0, or q = q t , is written from (33). With 17=0, 

* a 1 ~ n ° U 

tan 6= - 1 — — 
Its discussion is of minor interest. 

DRAG AND MOMENT 

Formulas for the pressure p all over the simple quadrics here treated are well known, for 
oblique as well as axial flow, and serve to find the drag and moment. For uniform flow the 
resultant drag is zero; its zonal parts can be found as heretofore. The moment about z is the 
surface integral of p(y dy dz — x dx dz), and generally is not zero. 

REGIONS OF EQUAL SPEED ABOUT OBLIQUE MODELS 

Compounding the velocities (31), (33) at any point in the ab plane, as in Figure 22, gives 
for q constant 

2 2 = [(l + m a )L7sin 0- (1 + m b ) V cos 0] 2 + [(1 -n a ) U cos 6+ (l-n b ) V sin 0] 2 = const. (41) 
In particular for q 2 = TP+ V 2 (41) gives 

(42) 



tan 0 = | (A± ^BC+A 2 )^ tan p. 

where K=V/U, and 

A = (l + m a ) fr+m^-(t^^(l-^ A 



n a (2-n a ) 



m b (2 + m b ). 




Pofar stream //he 
P' Stream axis 



Fig. 22.— Polar streamline and component velocities for uniform stream of inviscid liquid about oblique 

elliptic cylinder 



For an elliptic cylinder, as is well known. 



m a = 



b a+b 



which determines A,B, C, and thence 0 in terms of a' b'. Thus, for an endless elliptic cylinder 
of semiaxes a = 4, 6=1, yawed 10° to the stream, i. e., F/c7=tan 10° = . 1763, the graph of (42) 
has the form shown full line in Figure 23. This graph takes the dotted form when V= 0, q 0 = U. 



b a + b 



a a + b 
mb ~b'a' + b' J 



_a_ a + b 
Ub ~a' ar' + V 9 
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For a prolate spheriod of semiaxes a = 4, 6 = 1, yawed 10°, the graph of (42) is shown in 
Figure 17. 




Fhj. 23— Lines of steady flow, lines of constant speed and pressure, for infinite frictionless liquid 
streaming across endless elliptic cylinder. Dotted curve refers to stream parallel to x ; full-line 
curve q=q 0 refers to stream inclined 10° to x 



The two values of tan 0 in (42) are 

tan &=§(jt + y£(7+4 2 ), tan jBC+A% (43) 

from which are readily derived 

tan ^m^i^^f tan (ft + ft) = jffr . (44) 

(43) give the z-ward inclinations ft, ft, of the asymptotes of the curves g = g 0 . As can be 
proved, the inter asymptote angle ft- ft remains constant as K( = V/U) varies and the asymp- 
totes rotate through MKft + ft) about the c axis. 

Thus, with an elliptic cylinder, giving A, B, 0 their values at oo makes 

tan (ft tan (ft + ft) - ! ( 45 ) 

hence the asymptotes continue rectangular, as in Figure 23, while with varying angle of attack 
they rotate through J^(ft + ft). Or more generally one may show that da (ft - ft) = 0 .'. ft - ft = 
const. 

A similar treatment applies to the other figures of Table III. For all the cylinders the 
interasymptote angle is 90°; for the spheroids it is 2tan" 1 V2 = 109°-28' in the ab plane. 
Figure 17 is an example. If the flow past the spheroids is parallel to the he plane the inter- 
asymptote angle for the curves z = q 0 in that plane is obviously unaffected by stream direction. 
It is 90° for infinitely elongated spheroids; 109° -28' for all others. Excluded from the gen- 
eralizations of this paragraph are the infinitely thin figures, such as disks and rectangles edge- 
wise to the stream, that cause no disturbance of the flow. Passing to three dimensions, we 
note that the asymptotic lines form asymptotic cones having their vertex at the origin. 
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SUMMARY 

For an infinite inviscid liquid streaming uniformly, in any direction, past an ellipsoid or 
simple quadric: 

L The velocity potential at any confocal surface point equals the greatest tangential 
speed along that surface times the distance from the point to the surface's zero-potential plane. 

2. The tangential flow speed at said surface point equals the greatest tangential speed 
times the sine of the obliquity, or inclination of the local surface element to the equipotential 
plane. 

3. The normal speed at the point equals the greatest normal speed times the cosine of the 
obliquity. 

4. The locus of q = q 0 is a cup-shaped surface asymptoting a double cone with vertex at 
the center. 

5. The vertex angle of this cone is invariant with stream direction; for cylinders it is 
90°, for spheroids it is 2tan" 1 ^2 = 109° -28'. 

6. The velocity and pressure distribution are closely the same as for air of the same 
density, except in or near the region of disturbed flow. 

7. The zonal drag is upstream on the fore half; downstream on the rear half; zero on the 
whole. 

For the same stream, but with kinematic viscosity ?, if the dynamic scale is R = q 0 d/v, 
d being the model's diameter: 

8. The drag coefficient of a sphere is 24/R for R<.2; 28#- 85 + .48 for 0.2<#<200,000; 
and 0.5 for 10 4 <#<10 5 . 

9. The drag coefficient of an endless round cylinder fixed across stream is 8ir/R (2.002 — 
\og e R) for R<.5; approximately 9.4 R-*+1.2 for 0.5<#<200,000; 1.2 for 10 4 <i?<200,000. 

10. For 15,000<#<200,000 the drag coefficient of a round cylinder is 2.4 times that for 
a sphere. 
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TABLE III 

Flow functions for simple quadrics in stream V along y positive 
(For all shapes <f>= — (l + m b ) Vy, q t = — (l + m b ) V cos 0, q n =-(l — n b ) V sin 0) 



Shape 



Sphere 



2a'* 



Circular cylinder 



a' 2 



Elliptic cylinder 



Prolate spheroid 



Oblate spheroid 



a_ a + b 
b' a'+b' 



log. 



2c' 
l-e' 2 



1 — e 1 — e L 



COS co* 



e'Vl- 
" 1 + e 2 



Vi- 



p/2 — sin -1 *' 



— sin _1 c 



n b 



9l a +b 
a' a'+b' 



log ±±°- 2c' l - 2e ' 2 



-2e 2 



, 1 + e' 2 . _. , 

1+e 2 . ~ ' 
e / -sm J e 
VI — e 2 



* w is the longitude of any meridian; it is reckoned positive from y to z. 

TABLE IV 

Inertia factors k b for quadric surfaces in steady translation along axis b in Figure 20 



Elonga- 
tion E 


Effip. cyl. E=a/b 


Prol. spher. E=ajb 

l-e l-e 2 


Obi. spher. E=b/a 
jr. e— E sin l e 


**~ lag + « 2* 1 - 2 * 2 
1 — e 1 — e 2 




1. 00 


1. 00 


0. 500 


0. 500 


L 50 


L 50 


. 621 


. 384 


2. 00 


2. 00 


. 702 


. 310 


2. 50 


2. 50 


.763 


.260 


3. 00 


3. 00 


. 803 


. 223 


4. 00 


4. 00 


. 860 


. 174 


5. 00 


5. 00 


. 895 


. 140 


6. 00 


6. 00 


. 918 


. 121 


7. 00 


7. 00 


. 933 


. 105 


8. 00 


8. 00 


. 945 


. 092 


9. 00 


9. 00 


. 954 


. 084 


10. 00 


10. 00 


.960 


.075 


CO 


CO 


1. 000 


0 
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The numerical values in column 3 are given in Lamb's paper already cited; those in column 4 are given 
substantially by Doctor Bateman, Report No. 163 National Advisory Committee for Aeronautics, 1923. 
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Positive directions of axes and angles (forces and moments) are shown by arrows 



Axis 


Force 
(parallel 
to axis) 
symbol 


Moment about axis 


Angle 


Velocities 


Designation 


Sym- 
bol 


Designa- 
tion 


Sym- 
bol 


Positive 
direction 


Designa- 
tion 


Sym- 
bol 


Linear 
(compo- 
nent along 
axis) 


Angular 


Longitudinal 

Lateral 


X 
Y 
Z 


X 
Y 
Z 


rolling 


L 
M 

N 


Y >Z 


roll 


e 


II 
v 
w 


V 
Q 

T 


pitching 

yawing 


Z >X 


pitch 


Normal 


X > Y 


yaw_ 









Absolute coefficients of moment 

r _ L a -X-c -— 



Angle of set of control surface (relative to neu- 
tral position), 6. (Indicate surface by proper 
subscript.) 



4. PROPELLER SYMBOLS 



Pv, 



Diameter. 
Effective pitch 
Mean geometric pitch. 
Standard pitch. 
Zero thrust. 
Zero torque. 
p/D t Pitch ratio. 
V, Inflow velocity. 
V 9 , Slip stream velocity. 



T, Thrust. 
Q, Torque. 
P, Power. 

(If " coefficients" are introduced 
units used must be consistent.) 
v , Efficiency =TV/P. 
n, Revolutions per sec, r. p. s. 
N, Revolutions per minute., R. P. M. 

\27rr7i/ 



all 



<i>, Effective helix angle = tan 
5. NUMERICAL RELATIONS 



1 HP = 76.04 kg/m/sec.=550 lb./ft./sec. 
1 kg/m/sec. =0.01315 HP. 
1 mi./hr. = 0.44704 m/sec. 
1 m/sec. = 2.23693 mi./hr. 



1 lb. =0.4535924277 kg. 
1 kg = 2.2046224 lb. 
1 mi. = 1609.35 m = 5280 ft. 
1 m = 3.2808333 ft. 
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